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We investigate the possibility of obtaining chimera state solutions of the non-local Complex 
Ginzburg-Landau Equation in the strong coupling limit when it is important to retain amplitude 
variations. Our numerical studies reveal the existence of a variety of amplitude mediated chimera 
states (including stationary and non- stationary two cluster chimera states), that display intermit- 
tent emergence and decay of amplitude dips in their phase incoherent regions. The single-cluster 
chimera state is seen to co-exist with the spatially uniform state of the system whereas both types 
of two cluster chimera states exist in the parameter region where the uniform as well as the next 
higher traveling wave state are stable. These amplitude mediated chimera states may find useful 
applications in understanding spatio-temporal patterns found in fluid flow experiments and other 
strongly coupled systems. 
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Chimera states, spatio-temporal patterns of co-existing 
coherent and incoherent behaviour in an array of cou- 
pled identical oscillators, have received a great deal 
of attention in recent times [l|. First found by Ku- 
ramoto and Battoghtokh [2| from numerical investiga- 
tions of the weak coupling version of the non-local 
Complex Ginzburg-Landau Equation (NLCGLE), the 
chimera state has subsequently been studied for a va- 
riety of systems including two dimensional ones [3[ and 
those that have time-delayed coupling [4[ or those with a 
time-delayed feedback [5|. The phase only chimera state 
has been suggested as a useful paradigm to represent such 
curious phenomenon as uni-hemispheric sleep in certain 
mammals and birds, where during sleep one half of their 
brain is quiescent while the other half remains active |6| . 
Recently the phase only chimera states have also been ob- 
served experimentally in a chemical system [7[ as well as 
in an opto-electronic set up Q under controlled labora- 
tory settings. These past studies have however been con- 
fined to the weak coupling limit of the oscillator arrays 
where the amplitude variations have been ignored and 
only the dynamical behaviour of the phases have been 
considered. In many practical situations, such as in fluid 
flow representations, amplitude equations provide a more 
realistic description of the physical phenomena and have 
been widely employed to study the collective behaviour 
of such systems. An interesting question to ask is there- 
fore whether spatio-temporal patterns corresponding to 
chimera states can exist for the strong coupling limit. In 
this paper we have addressed this question and we report 
the numerical discovery of chimera solutions for the NL- 
CGLE equation in the regime where amplitude effects 
matter. In contrast to the classical chimera states found 
for the phase only systems, the present ones display am- 



plitude activity in the incoherent region of the solution 
in the form of intermittent emergence and decay of am- 
plitude dips and we classify them as amplitude mediated 
chimeras (AMCs). The phases of the oscillators in the 
incoherent region continue to have a random distribution. 
These states bear a close resemblance to the simultaneous 
appearance of laminar and turbulent regions in Couette 
flow studies @ and may have wider applications to other 
strongly coupled systems. 

Our model system is the well known one dimensional 
NLCGLE |1(| that has been extensively studied in the 
past in the context of applications to various physical, 
chemical and biological phenomena |ll| . 
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where W{x,t) = A(x,t)exp[i(j)(x,t)] is a complex field 
quantity with A(x, t) and 0(x, t) representing the am- 
plitude and phase respectively. C\, C 2 and K are real 
constants characterizing the linear and the non-linear 
dispersion and the coupling strength respectively. The 
non-local mean field W(x, t) is given by 



W(x,t) 



G(x-x')W{x-x',t)dx'. 
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where the coupling kernel G(x — x') is taken to be of an 
exponentially decaying nature and its normalized form is 
given by, 
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Here k > is the inverse of the coupling range and is a 
measure of the nonlocality of the coupling. The space co- 
ordinate is made dimensionless by normalizing it with the 
system length L and hence the system size extends from 
— 1 to 1 . Very large to very small k in comparison to the 




FIG. 1: (Color online) The upper panels exhibit snapshots of the stationary state spatial profiles of the phase (<f>) and the 
amplitude \W\ (multiplied by 10 and shown in blue) for (a) a single cluster AMC with K — 0.40, Ci = —0.6, Ci — 2.0; (6) 
a 2-cluster AMC with K = 0.40, Ci = -4.0, C 2 = 0.5; and (c) a phase-only chimera with K - 0.05, Ci = -0.9, C 2 = 1.0 
respectively, k is fixed at 2 for all the cases. The lower panels (d — /) show the corresponding spatial profiles of the long-time 
average of the order parameter amplitudes (R) and the frequencies (w =< <j> >)• The computed value of A is marked with a 
horizontal dashed line (in black). 



system size (which is 2) steers the system through local- 
to-intermediatc-to-global coupling regimes. The magni- 
tude of the constant K determines the coupling strength 
and for K « 1, Eq.{T]) can be reduced to a non-local 
evolution equation for the phase function (j)(x, t) that has 
been the subject of several past studies for the classical 
chimera solutions. We have carried out extensive nu- 
merical explorations to seek chimera solutions of Eq.([T]) 
and have discovered a variety of such states over a wide 
range of parameter space. Broadly they consist of sta- 
tionary one cluster and two cluster chimera states and 
also a non-stationary {breather) variety of the two clus- 
ter state. The two coherent regions of the two cluster 
chimera states have opposite phases and are separated by 
incoherent regions. Multi-cluster phase coherent regions 
have also been observed before in time delayed systems J4J 
and more recently by Zhu et al. [12| in the weak coupling 
limit of the NLCGLE. The major difference of the present 
solutions from their counter parts of the phase only sys- 
tems is that these amplitude mediated chimera states 
have significant temporal variations of the amplitude in 
the incoherent spatial regions. These regions show in- 
termittent emergence and decay of amplitude dips which 
in some cases can resemble amplitude hole (defect) so- 
lutions. In looking for chimera states, our choice of the 
system parameters have been guided by earlier investi- 
gations of the NLCGLE including in the weak coupling 
limit. Thus we have chosen two values of K, namely 0.05 
and 0.4 to represent weak and strong coupling cases re- 
spectively. The value of k has been chosen to be equal 
to 2 so that kL = 4 which is the same as chosen by 



Kuramoto and BattogtokQ (who had k — 4 and L = 1) 
where L is the system size. The values of G\ and C^ have 
been varied over a wide range. In FigfT]we show a typi- 
cal snapshot of the amplitude mediated one cluster and 
two cluster solutions in panels (a) and (b) respectively 
whereas in panel (c) for comparison we display a classical 
phase only chimera obtained in this case from Eq.(TT|) in 
the weak coupling limit by taking a very low value of K. 
Notice that for the classical chimera state the amplitude 
fluctuations are negligible, justifying their neglect in the 
weakly coupled limit of the NLCGLE. For the classical 
(phase only) chimera the value of C\ is —0.9 and G<x = 1 
corresponds to a = tan _1 (C 2 - Ci)/(1 + CiC 2 ) = 1.52 
where a is the phase-shift parameter in the weak cou- 
pling limit (3|. For the AMCs the amplitude variations 
in the incoherent region are quite significant, sometimes 
dipping close to zero values corresponding to traveling 
hole like solutions. Our simulations have been done with 



the XPPAUT 13| package with 201 discrete oscillators 
equally spaced on a ring. We have carefully checked our 
numerical results to rule out finite size effects. The na- 
ture of the AMC is found not to change when, for exam- 
ple, we change the number of oscillators from 201 to 401. 
The characteristics of the AMC remain the same and it 
does not exhibit any transient nature or tendency to col- 
lapse. The initial conditions of our simulations consist of 
slightly perturbed uniformly spaced phases from to 2tt 
with unit amplitude. 

The stationary patterns of these AMCs can also be un- 
derstood in terms of a complex order parameter defined 
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(4) 
where i?(x, t) is the amplitude of the order parameter, 
Q(x, t) is the mean phase and 8 = <fi + fit is the relative 
phase defined in the frame rotating with the angular fre- 
quency Q and amplitude A of the coherent segment of 
the chimera. Using Eq.Q, and separating Eq.([T|) into 
its real and imaginary parts, one can get, 
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By restricting to time stationary solutions, we get: 
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The absolute value of the right hand side of Eq. ^ cannot 
be greater than 1 and this puts a condition on the mag- 
nitude R of the order parameter, namely (R(x) > |A|) 
in any coherent segment in space where 
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We obtain the amplitude A and the frequency $7 of the 
coherent segment from the simulations and compute the 
value of A using Eq.©. In FigfTJd)-(f) we have plot- 
ted the time averaged profiles of R(x) for the chimera 
states corresponding to the snapshots in Figflja)-(c) re- 
spectively. The horizontal line in each figure marks the 
computed A value. As can be seen the results are in 
good agreement in that the coherent segments found in 
FigfT] correspond to regions where the condition R > A 
is satisfied. As a further check on the nature of the col- 
lective state we have also plotted the average frequency 
profiles in FigQ] (d)-(f) which all show the typical sig- 
nature of chimera states, namely a constant frequency 
in the coherent region (flat profile) and a peaked pro- 
file in the incoherent region [2J. For the non-stationary 
breather state the order parameter is no longer a constant 
quantity but shows a periodic temporal variation. This 
is shown in Figl^k for the 2-cluster AMC state. Fig(5b 
shows the spatio-temporal pattern of the phases 4> after 
the transients. We notice the oscillating nature of the 
the two cluster chimera. To locate the existence regions 
of the AMCs we have next carried out a linear stabil- 
ity analysis of plane wave solutions of Eq.([T]), that are of 
the form, W%(x,t) = ake K*k*-u k t) and that sat i s f y tne 

following dispersion relations, 
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FIG. 2: (Color online) (a) The temporal patterns of the am- 
plitude \W\ of an oscillator and the modulus R (the lower 
curve in blue) of the complex order parameter, (b) The spatio- 
temporal pattern of the phases with time running downwards 
from T — 1800 to 2000. The parameter values for this simu- 
lation are k = 2, K = 0.4, Ci = -8 and C 2 = 0.95. 



where K' = K(l - I K , k ) and I K>k = j\ G(x')e™ kx ' dx' . 
Perturbing these equilibrium solutions by writing W k = 
[1 + u(x, t)]W%, where u(x, t) = J2 n u n (x, t)e xt e l7Tnx and 
substituting in the linearized form of Eq.([T]), we can get 
a variational equation for u n , 



du n 
~dt 



[1 + iuo k - 2(1 + iC 2 )a\ + K{1 + id) 

x {I K , n+k - 1)K - (1 + iC 2 )a\u n (10) 



Eq. (fT0|) and its complex conjugate yield a 2 x 2 matrix 
M , whose eigenvalues are determined from the following 
quadratic characteristic equation |14j . 

\M - XI\ = A 2 + (n + ir 2 ) A + Pl + ip 2 = (11) 



where / is the 2x2 unit matrix and ?*i = — (a + e), r 2 = 
— {b + f),pi = —bf + ae — c 2 — d 2 and p 2 = af + be, which 
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FIG. 3: (a)Stability diagram for the synchronous state (k — 0, solid curve) and for k — 1 traveling wave (TW) state (dashed 
curve) in G\ — K plane with k = 2 and C2 = 1. The region below the solid curve denoted by U is unstable to both the states 
whereas the region below the dashed curve not lying within U is stable to k — 1 TW state. S denotes the stability region for 
the synchronous state and B marks the bistable region where both the synchronous and the TW states are stable. The phase 
approximation is valid in the vicinity of the dotted horizontal line at K — 0. (b) and (c) The stability diagrams similar to 
panel a but in the phase plane of C\ — C%. The coupling constant K is chosen to be 0.05 for the panel b and 0.40 for panel c. 
The points denote the parameter values for which various chimera states were obtained through simulations. The simulation 
results for the four filled plot-markers (one in panel b and three in panel c) are shown in this paper. The circular plot-markers 
correspond to phase chimera in panel b and the single cluster AMCs in panel c and the square ones denote the 2-cluster AMCs. 
The extreme left diamond plot-marker in panel c corresponds to a 2-cluster breather AMC. 



in turn are expressed in terms of system parameters as: 
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Eq. ljlip determines the eigenvalues A for a perturbation 
with a wave number n. The real part of the eigenvalues of 
Eq. tfTTj) should all be negative for all integer values of n. 
This gives us an important quantity determining the sta- 
bility of plane wave oscillations against the perturbation 
of a given wave number n. 
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Positive $ implies stability and negative $ instability. 
From our numerical analysis we find that the lowest mode 
number perturbation [n = 1) is the first one to get desta- 
bilized and therefore determines the marginal stability 
curve. We fix n = 1 for our further stability analysis. 
For the uniform (k = 0) state, we are able to get a sim- 
ple analytic form for the marginal stability curve, namely, 



$(*, 0, 1) = 1 + CyC 2 + K{1 + Cl2) (1 - 7) (14) 



where the term 7 = I K 1 
non-locality in the system. 



T , ■? accounts for the 



The stability condition is $ > and the expression for 
$ given by Eq. (|14[) reduces to the well known Benjamin- 
Feir- Newell criterion \ + C\C 2 > for 7 — > 1 correspond- 
ing to large k i.e. local or diffusive coupling. In the global 
limit (7—^0 for small k ), the stability condition reduces 
to that reported in earlier works [15(. FigO^. shows the 
stability diagram for the uniform state (fc = 0) as well as 
for the fc = 1 traveling wave state in the C\ — K phase 
space where C 2 is fixed at 1 and k at 2. The phase only 
models are valid near the dotted line at K — 0. Sim- 
ilarly Fig(3b,c show the stability diagrams in C\ — C 2 
phase space for two different values of K but the same 
value of k. The location of the chimera states are marked 
by different point symbols (filled circles, squares and di- 
amond) on these stability diagrams. We draw attention 
to the regions S and B of the stability diagrams where 
chimera states are found. We note that the single clus- 
ter chimera points are located where the uniform state is 
stable whereas the two cluster chimera points lie in the 
region where the uniform as well as the k = 1 traveling 
wave states are both stable. This bi-stability property of 
the NL CGLE leading to the co-existence of the uniform 
state and the single cluster chimera state has been ob- 
served for the classical chimera state as well and appears 
to be preserved in the strong coupling limit. Likewise the 
association of multi-cluster chimeras with correspond- 
ing wave number traveling waves (as noted previously 
in phase only systems [J, |l6( ) is once again noticeable for 
AMCs. 

In conclusion, wc have studied the NLCGLE system 
in the strong coupling limit and found a new class of 
chimera states where the incoherent regions display sig- 



nificant amplitude fluctuations. These amplitude medi- 
ated chimeras can be of the stationary kind (with a single 
or two cluster configuration of coherent regions) or have 
an oscillatory nature. Our detailed numerical investiga- 
tion have also marked out the existence regions of these 
hybrid states in the reference frame of the stability dia- 
gram of the uniform state and the k = 1 traveling wave 
state of the NLCGLE. These states not only complement 
the previously found phase only chimera states but also 
extend the applicability domain of such hybrid states to 
physical systems that are better represented by full blown 
amplitude equations such as the NL CGLE. Some systems 
that come to mind in this context are fluid flow simula- 



tions/experiments where the simultaneous appearance of 
laminar and turbulent regions have been observed |9[ and 
neuronal networks displaying bump states where a subset 
of neurons fire in synchrony while others fire incoherently 
[17| . The discovery of these novel states also opens up a 
number of interesting future areas of investigation includ- 
ing a study of their stability, delineating their linkages to 
other coherent solutions of the NLCGLE such as trav- 
eling waves and holes and exploring their existence for 
other forms of the coupling kernel. 
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